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A number of Ce and U based heavy electron metals have recently been discovered to show "non-Fermi liquid 
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We study a new exchange interaction in which the conduction electrons with pseudo spin Sc = 
3/2 interact with the impurity spin Si = 1/2. Due to the overscreening of the impurity spin by 
higher conduction electron spin, a new non-trivial intermediate coupling strength fixed point is 
realized. Using the numerical renormalization group (NRG), we show that the low-energy spectra 
' are described by a non-Fermi liquid excitation spectrum. A conformal field theory analysis is 

^\ I compared with NRG results and excellent agreement is obtained. Using the double fusion rule to 

0^ ■ generate the operator spectrum with the conformal theory, we find that the specific heat coefficient 

and magnetic susceptibility will diverge as T~^^^, that the scaling dimension of an applied magnetic 
f— I , field is 5/6, and that exchange anisotropy is always relevant. We discuss the possible relevance of 

• our work to two-level system Kondo materials and dilute cerium alloys, and we point out a paradox 

in understanding the Bethe-Ansatz solutions to the multichannel Kondo model. 
QQ , * Present address. 
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■ physics," in such anomalous properties as: i) a specific heat coefficient linear in logT, T the temperature; ii) resistivity 
I which approaches its T = value as with a close to 1; iii) a magnetic susceptibility x(^) which cither diverges 
as a logarithm or weak power or saturates as x{T) — x(0)[l — A\/T] In the case of UBcia and CeCu2Si2, a 

[ superconducting instability arises within this exotic normal state . It is an open question whether these anomalous 
properties arise from distortions of Fermi liquid physics such as arise in the vicinity of a quantum critical point 0| 
or due to a disordered distribution of Kondo scales [g| , or whether the anomalies signal the appearance of an entirely 
new fixed point regulating the low temperature physics of these materials. 

In this latter regard, it is prudent to study impurity models, such as the overcompensated multi-channel Kondo 
Q ' model 1^ which have very plausible relevance to Ce and U based systems In this model, M identical spin 1/2 

conduction bands exchange couple to a single spin 5*/ impurity with the condition M/2 > Si, so that there is more 
conduction spin than needed to fully compensate the impurity, which leads to an instability of the strong coupling 
(exchange coupling J — > 00) fixed point, as indicated in Fig. 1. Because the weak coupling (J = 0) fixed point is 
^ unstable always due the Kondo effect, this implies the presence of a non-trivial fixed point at intermediate coupling. 
This non-trivial fixed point corresponds to a degenerate ground state and non-Fermi liquid excitation spectrum-no 
1:1 map exists of the states to a non-interacting fermion system. In particular, the charge and spin relations relevant 
for free fermions do not apply for these systems where spin/charge/channcl separation occurs j^. For more than one 
impurity, this fixed point cannot remain stable at low temperatures |]l0|] , due to the overlap of the polarization clouds 
around each site which have diverging length scales as T —^ 0. It is conceivable that the impurity models may regulate 
the physics of alloys with finite concentrations of Ce and U ions or even concentrated compounds over a finite range of 
temperatures before ultimately giving way to new fixed points as T ^ 0. This picture is supported by recent results 
in infinite spatial dimensions |lj[,|l^ . 

Recently, impurity models containing a single high spin electron band have been derived which have been argued 
to show non-trivial physics. These models have been argued to describe various unstable fixed points of the two-level 
system Kondo model (describing electron assisted tunneling of atoms in a double well) |l^], and a Ce impurity in a 
metal when finite Coulomb interaction is retained ||l^,|l5| . The simplest example of this kind of model is for impurity 
spin Si — 1/2 and conduction spin Sc = 3/2. Although only one channel of conduction electron is present, this model 
gives a nontrivial fixed point for the same reasons as the multi-channel model, as illustrated in Fig. 2. The point is 
that a single large spin Sc\tot = 2 is coupled to the impurity as the lowest energy state in strong coupling. Again, 
strong coupling is unstable, and weak coupling is unstable, so that a non-trivial intermediate coupling fixed point 
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must exist. This kind of fixed point has been studied with a perturbative scaHng analysis |13|-[L5[, and more recently 
with conformal field theory analysis |pj|] . 

In this paper we present a comprehensive numerical renormalization group analysis of the Sj = 1/2, Sc = 3/2 
single channel Kondo model combined with a conformal field theory study of the energy and operator spectrum. We 
confirm the existence of a non-trivial fixed point with a non- Fermi liquid excitation spectrum. We find good agreement 
between finite size spectra predicted by the conformal theory with those obtained from the NRG, both from the group 
theoretic basis we use for our conformal theory and that of Ref. ||l6| which uses a different basis. Using the double 
fusion rule of Affleck and Ludwig , we work out the spectrum of primary field operators for the model, and find that 
the magnetic susceptibility and specific heat coefficient in this model will diverge as T~^/^ for T — > 0, that the scaling 
dimension of the applied magnetic field is 5/6, and that exchange anisotropy will always be relevant in this model. In 
concluding, we briefly discuss the possible relevance of this model to intermediate temperature regimes of two-level 
system Kondo materials and dilute cerium alloys. Finally, we discuss a separate point of theoretical interest concerning 
a paradox about the connection of this model to the Bethe-Ansatz solutions of the overcompensated multi-channel 
Kondo model ^^j^. 

Our paper is organized as follows: In Section II we discuss the model Hamiltonian and the numerical renormal- 
ization group methodology used to study it. We also discuss the results from the numerical renormalization group 
calculations. In Section III we develop the conformal field theory concepts required to study this model and calcu- 
late the finite size spectra and primary operator spectrum, and we compare the conformal theory results with the 
numerical renormalization group results. In Section IV we summarize and conclude. 



II. NUMERICAL RENORMALIZATION GROUP 



The numerical renormalization group (NRG) technique was invented by K. G. Wilson to study the Kondo problems. 
This method defines the renormalization group (RG) Hamiltonians and diagonalizes the RG Hamiltonian numerically. 
The essence of this method lies in the logarithmic discretization of the conduction band and in transforming the 
conduction electron Hamiltonian into a tridiagonal form. Hence we can visualize the NRG Hamiltonian as a linear 
array of sites with the hopping integrals between sites. The impurity lies at the = — 1 site and the conduction 
electrons sit at the sites A^ = 0,1,2,- • •. 

Since this method depends on direct numerical diagonalization, the numerical work requires a huge memory storage 
compared with, e.g., the NCA. The consideration of all the possible symmetry of the Hamiltonian is essential in 
reducing the computer memory. The NRG is defined by a set of energy eigenvalues and the reduced matrix elements 
of the electron operators between the energy eigenstates. 

The NRG approach was the first to show the crossover physics of the Kondo effect from the high temperature regime 
to the low temperature fixed point. K. G. Wilson solved the simple Kondo exchange interaction model, calculated 
the magnetic susceptibility, and showed the RG flow from the high temperature fixed point to the low temperature 
fixed point. Later, the NRG was applied to the simple one-channel Anderson model with the conduction electron 
spins S'c = 1/2 interacting with the impurity system of 5*/ = 1/2 |^^. The two-channel Kondo model was also studied 
using the NRG |2^. The NRG method was able to show that the two-channel Sc = Sj = 1/2 Kondo model leads 
to the non- Fermi liquid ground state with a non-trivial finite exchange coupling at the fixed point. The one-channel 
two-impurity and two-channel two-impurity Kondo problems have been also studied using the NRG [^2yiC| ] and the 
low energy levels from the NRG have been compared to the finite energy spectrum results of conformal field theory 
to verify the fusion rule hypothesis used in the CFT approach . 

The NRG is a unique approach for studying stability of various fixed points of the Kondo problems. In contrast to 
other techniques, the NRG can provide the RG flow diagrams. In an earlier work, stability of the two-channel non- 
Fermi liquid fixed points were studied with the NRG against the channel anisotropy, the exchange coupling anisotropy, 
and the external magnetic field . This method was extended to the studies of the thermodynamics ||2^ and the 
dynamical properties [p5|-[28| like the transport coefficients and the local spectral functions. 



A. The conduction electrons with Sc = 3/2. 

Here we apply the NRG method to the case that the conduction electrons with Sc = 3/2 interact with the impurity 
spin Si = 1/2. Derivation of this exchange interaction is presented in Ref |l5|. Higher conduction electron pseudo 
spin can overcompensate the impurity spin leading to the non- Fermi liquid fixed point. In this section, we are going 
to show that this model leads to new non- Fermi liquid fixed point using the NRG. 
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1. Logarithmic Discretization. 



The Model Hamiltonian reads 



H = H, 



cb 



Hi, 



E 



Hi=J SciO) ■ Si. 



(1) 
(2) 

(3) 



Here the conduction electron spin Sc = 3/2 and the impurity spin Sj — 1/2. We will assume the flat square DOS 
for the conduction band. Following the standard transformation using the logarithmic discretization the NRG 
Hamiltonian reads 



oo 



1 - A-"-i 



in 
foa 



v/(l-A-2«-i)(l-A-2«-3)' 



de Cf 



%/2 J-i 



(4) 
(5) 

(6) 

(7) 



Here L is the canonical matrix representation of spin J = 3/2. This Hamiltonian can be envisaged as a linear array of 
sites with an impurity sitting at the site n = —\. Only the electrons at the origin (n = 0) interact with the impurity 
spin. This tight-binding form of Hamiltonian has a nearest neighbor hopping integral whose magnitude depends on 
the distance from the impurity site and is progressively reduced as the sites move away from the origin. 
We may define the Hamiltonian up to iV-th site by Hn 



1 



JV-l 



A Cri [ fnafn+la + fl+iafna ] + Hq, 

Ho - 2N{0)J foo.Lo.pfo0 ■ S. 



n=0 



(8) 
(9) 



We can diagonalizc this Hamiltonian step by step. The strategy is very simple. When the Hamiltonian upto the A^-th 
site is assumed diagonalized, we add the N + 1-th electron site and diagonalizc the extended Hamiltonian again. We 
continue this process until the structure of low lying energy levels does not change, or, when the stable fixed points 
are accessed. Since the hopping integral decreases exponentially with site number TV, rapid saturation is expected for 
a reasonable value of A. To make the lowing lying energy of order one, we redefine the Hamiltonian, 



Kn 



1 + A 



— A(^-i)/2 Hn 



^ A-"/2 e„ [ + fi+iJno. ] + 2iV(0) J fl^L^pfop ■ S 



,n=0 



Kn+1 = A^/^ Kn + E [ f Naf N+la + flj+iafNa ]■ 

a. 

The last equation defines the renormalization group of the Kondo Hamiltonian. 

2A-1/2 



Ko = Z — i?0 = ^ foa^apfop ■ S, 

aP 



1 + A 
4A-1/2 
1 + A^ 



- JiV(O). 



(10) 

(11) 

(12) 
(13) 
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2. Symmetries. 



Now we show that the axial charge and the total spin operators commute with the Hamiltonian. We note that the 
axial charge may be more commonly called the isospin, and that this nomenclature is peculiar to the NRG literature. 

The conduction band Hamiltonian is invariant under the unitary transformation /„„ — > Ua/sfn/d- Including the 
exchange interaction, the total spin is conserved 

J^Y.fnJo.0fnl3+S. (14) 
n 

It can be easily shown that the total spin of the conduction electrons and the impurity spin commutes with the 
conduction band Hamiltonian 

[-ffcb, /]=0. (15) 

Also this spin operator commutes with the usual Kondo exchange interaction. 
The Hamiltonian is invariant under the unitary transformation: Fna ^ UFna 

fna 



This invariance leads to the definition of axial charge symmetry operators 

n 

Q- = [Q+]t ^ fc^(-l)" ec^p /„/5/„«, (18) 



f'! f _ i 

JnaJna ^ 



(19) 



Here eap is an antisymmetric tensor to be determined. These axial charge operators commute with the conduction 
band Hamiltonian. 

[ ifcb, QM - 0. 

We will determine the prefactor k and the antisymmetric tensor by demanding that these axial charge operators 
satisfy the angular momentum commutation relations. First they satisfy 

[Q^Q±]=±g±. (20) 

We demand 

[Q+,Q-]=2g^ (21) 
Then we find the following conditions to be satisfied 

= [^Xp^-Sap- (22) 



Hence we may choose, e.g., for Sc — 3/2, 



fc = i; e^a'^'^iay. (23) 



Thus 

oo 
n=0 

The above axial charge and angular momentum operators commute with each other and with the ordinary exchange 
interaction. 
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We can also generalize the axial charge operators for the conduction electron band with any half-integer spin. The 
angular momentum operators are self-evident. 



(25) 



Q- = 



fnafn 



2 ^ 



{Nn- 



(26) 
(27) 



Thus the angular momentum operators and the axial charge operators commute with both the conduction band 
Hamiltonian and the ordinary Kondo exchange interaction. They also commute with each other, and hence provide 
good quantum numbers for our numerical calculations. 



3. Irreducible symmetry tensor operators. 

In this section, we will construct irreducible tensor operators at each site which are building blocks for new symmetry 

states in the NRG work. 

The angular momentum operators are 

= E " fLfNa, (28) 
a 

= V3 /lf3/2/jVl/2 + 2 flfi/2fN-l/2 + fl,_i/2fN-3/2, (29) 

= \/3 7^1/2/^3/2 + 2 /^_i/2/jVl/2 + flj_^/2fN-l/2- (30) 

And the axial charge operators are 

^E (31) 

Q. 

Qn = ( /iV3/2/jV-3/2 ~ /Ari/2/iv-l/2 )' (^2) 

QN = [Qtr (33) 

The electron operators fj^^ and the hole operators h]^^ = (— 1)-'^+3/2+« Jj^- are irreducible tensors of rank 3/2 in 
the angular momentum sector. That is, 

[ /L ] = « /L, (34) 

[ /L ] = /L+i = V-K.'^ + 1) - «(« + 1) /L+i. (35) 

[ S],, /L ] = i.-i,. /L-i = Vsis + l)-a{a-l) fj,^_„ (36) 

[ S},, /^L ] = « h^Nc^ (37) 

[ h^No. ] - ^L+i = V5(s + l)-a(a + l) (38) 

[ h^Na ] = ia-i,a '''jVa-i = \/ s{s + I) - a{a - 1) h^j^^_^. (39) 

Here i+ and L~ are the standard S = 3/2 angular momentum matrices. Hence we can apply the Wigner-Eckart 
theorem when we calculate the matrix elements between two different angular momentum operator symmetry states. 
In the axial charge case, the electron operators and the hole operators satisfy the following commutation relations 

[Qn, /iL] = -^/ik' (40) 

[ Qn, fL] = l fL, (41) 

[Qtf, -/^k]=/L. (42) 

[ Ql,, fL ] = 0> (43) 

[Qn, fL]=-hL, (44) 

[ Qn, ] = 0. (45) 
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Thus wc can sec that the pair {flf^j ^^]va) form a tensor of rank 1/2 in the axial charge sector. Thus we can apply 
the Wigncr-Eckart theorem to these conjugate pair operators. 

Wc now construct the irreducible symmetry operators, Tn{QR; JM), at the site A'' using the above operators. At 
each NRG site, the following three types of irreducible symmetry tensor operators can be developed and they generate 
16 new symmetry states: 

1. Q = 1, J = 0: We choose the vacuum state and operate the axial charge raising operator to find all states 

Tw(li;0) = l (46) 
Tjv(10;0) = 4s Q+ Tjv(ll;0) 



^/2 

(-IT 
V2 



(/jv3/2/jv-3/2 ^^1/2/^-1/2) 
T;v(ll; 0) = ^ 0^ T^(10; 0) = -/L/2/m/2/Ll/2/L3/2- (48) 

2. Q = 1/2, J = 3/2: First we define the Q = 1/2, R = —1/2 states and find the R = 1/2 states operating the 
axial charge raising operator on them 

I^iv(^(-^);^m)=/^„. (49) 

We already showed that these operators satisfy the proper commutation relations with the angular momentum 
operators. The R = 1/2 operators are found by operating the axial charge raising operator on the above 
operators 

Tn{\\: Im) = Q+ Tj,{\{-\)- ^m) = h)^^ T^{11;Q), (50) 
/^k = (-1)^+^/^+™ fN^. (51) 

As noted before, the hole operators h}^^^ satisfy the tensor (of rank 3/2) commutation relations with the angular 
momentum operators. 

3. Q = 0, J = 2: The irreducible operators can be found using the commutation relations 

T^(0; 22)^/^^3/2/^1/2, (52) 
T^(0; 21) = \[ S],, r;v(0; 22) ] = /^3/2/Li/2' (53) 
Tjv(0;20) = -^ [S],, T;v(0;21)] 

(/jV3/2/jV-3/2 + /jVl/2/lr-l/2) ' (5^) 



V2 

Tn{0; 21) = ^ [ S^, T^iO; 20) ] = /l^i/a/Ls/s' (55) 
T^iO; 22) = i [ 5^, T^iO; 21) ] = /Li/2/L3/2- (56) 



4- Construction of the NRG Hamiltonian and the reduced matrix elements of the discretized conduction electron operator 

16 new irreducible symmetry electron operators are generated at each NRG site. These new operators belong to 
one of three tensor operators Tn{QR; JM) with rank Q in the axial charge sector and with rank J in the angular 
momentum sector 

Tn{00;2M), Tn{^R;^M), Tm{1R;00). (57) 
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These irreducible tensor operators are basic building blocks for the construction of new symmetry sector. 

In constructing the NRG Hamiltonian in each symmetry sector or in calculating the reduced matrix elements of 
the discretized conduction electron operators, we have to calculate the following matrix elements 



( T^(li?i;00) 
T*(00;2Mi) 
, T^(ii?i; |Mi) , 



rjv(ii?2;00) 

fL { TNm2M2) 

rjv(ii?2; IM2) . 



According to the Wigner-Eckart theorem, the non-vanishing combinations are 



rjv(ii?2;00) 

r^(00;2Af2) 



r^(li?i;00) 



/krjv(^i?2;^M2). 



r^(00;2Mi 

Guided by the Wigner-Eckart theorem, the explicit calculation gives 



I <^^il^^'1^2> (Ml, a) 
,111 



r^(ii?i;^Mi) fL T^(00;2M2) = (-1)^+^ ^ < ^Mi|^a;2M2 >, 

T* (00;2Afi) fL Tn{\r2; ^M^) = V2 {R2.-\) < 2Ah\^a; h'h > . 



(58) 



(59) 



(60) 
(61) 

(62) 
(63) 



Here {x,y) is the shorthand notation for the Kronecker delta function S^ y. Using the above equations and the 
Wigner-Eckart theorem, we can construct the NRG Hamiltonian matrix elements at each symmetry sector which are 
expressed in terms of the reduced electron opreator matrix and the 6 ~ j symbols. 

As mentioned before, the numerical renormalization group is reduced to book-keeping of two sets of energy eigen- 
values and reduced matrix elements in the Kondo problems. The reduced matrix elements are required for the 
construction of the NRG Hamiltonian in the next NRG site. Hence it is essential to keep track of two data sets. In 
the next section, we are going to present our NRG results, and give a physical interpretation. 



B. NRG results 



In this section, we are going to present the results of NRG calculation for the exchange interaction of the conduction 
electrons with spin Sc = 3/2 with the impurity spin Sj = 1/2. 

First we show the low energy spectra in the weak coupling limit with the initial coupling g — 10^^° (free electron 
limit). We used the discretization parameter A = 5 and the truncation energy 10. In this extremely weak coupling 
limit, the low lying excitation energy spectrum should be described by the free electrons. As shown in Fig. ^ the low 
lying energy levels are evenly spaced. 

When we increase the bare exchange coupling, we should be able to observe the crossover from the weak coupling 
fixed point to a nontrivial non-Fermi liquid fixed point. In Fig. ^ we display the RG flow with initial exchange 
coupling g = 0.01. We can observe the clear crossover from the weak coupling fixed point to new fixed point. The 
low lying excitation energy spectrum at this new fixed point is completely different from that of free electron at the 
weak coupling fixed point. That is, this new fixed point is non- Fermi liquid fixed point. ^From the Fig. ^ we can see 
that degeneracies at weak coupling fixed point are lifted after crossover at new non- Fermi liquid fixed point. 

Now we present the three sets of figures with initial couplings ranging from weak coupling (see Figs. |^), coupling 
close to the fixed point(5* ~ 0.2 form the third-order scahng analysis, see Figs. ||), and strong coupling(see Figs. ^). In 
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these cases, we used the same model parameters as above. The low energy structure is completely different from that 
in the zero coupling limit (free electron). That is, this interaction model generates new non-Fermi liquid fixed point. 
Three different initial couplings flow to the same fixed point as can be seen from the fixed point low energy spectra. 
The ground state is alternating between {q = 0, j — 1/2) sector and {q = 0,j = 3/2) sector with the conduction 
electron screening shell (NRG site N). For odd N, the {q = 0, j = 1/2) sector is the ground state and vice versa. 
Here {q — axial charge, j = angular momentum). 

In summary, we have shown using the NRG that the one-channel Sc — 3/2,5'/ — 1/2 exchange interaction model 
generates non-trivial fixed point which is completely different from the Fermi liquid fixed point. 



III. CONFORMAL FIELD THEORY ANALYSIS: LOW-ENERGY SPECTRUM 

In this section, we carry out a conformal field calculation (3^,^ determining the first few states and eigenvalues 
in the spectrum of the low-energy fixed-point. In contrast with the NRG procedure, the GFT approach starts out 
with the model Hamiltonian expressed in real-space coordinates. Following Affleck and Ludwig we introduce 
left-moving field operators ^pa(x) {a = ±3/2, ±1/2) and rewrite Eqs. (2) and (3) as 

/oo 
dx i>i{x)dMx), (64) 
-oo 

and 

Hi = VFpJMOyL^pMO) ■ Si, (65) 

where p — 1/D is the density of states. 



A. Free electrons 



We must next find the Sugawara forms |^ equivalent to Eqs. (|^ and (|65|). These are harmonic oscillator-like 
Hamiltonians quadratic in certain "current" operators of the system whose spectrum can be determined by the relevant 
ladder-operator algebra. We follow the logic of Affleck and Ludwig first, decompose the free electron Hamiltonian 
in a Sugawara form which separates out the spin currents from other currents (e.g., charge) of the system. Then 
couple the electron spin to the impurity and generate the new spectra. The interaction with the impurity spin is 
"absorbed" into the spin sector of the Sugawara form Hamiltonian through a simple completing of the square. The 
excitations of the interacting system are found from those of the non-interacting system through application of the 
"fusion rule" |^] which, loosely, corresponds to simple addition of the impurity spin to the conduction spins subject 
to the constraints of the Kac-Moody spin algebra. 

We first consider the case of anti-periodic boundary conditions for the conduction states which corresponds to the 
odd-TV NRG iterations, and gives rise to a non-degenerate free electron ground state. The Sugawara Hamiltonian is 
constructed from the Fourier components of the conserved currents — axial charge (or isospin) and spin. We therefore 
consider a segment of length 2£, with anti-periodic boundary conditions, and define a Fourier sequence for each 
component of the axial charge: 



exp(— i7rna;/£) dx, 



(66) 



and 



Qn^{Qn) V'g/aWV'la/ala;) - V'I/2(a;)'/'li/2(a;) exp(-i7rna;/£) dx. 

Expressed as combinations of the conduction band operators c^, these components are 



E 



(67) 



(68) 



and 
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Qn ^ {Qn) — ^ (ce,3/2cL£_„TrDF/2£,-3/2 ^l,l/2'^^-€-mrVF/2e,-l/2) ' (^^) 
e 

Likewise, for the spin current, we have 

J„ = / ^pl{x)Lat3i^i3{x)exp{-iTTnx/i),dx, (70) 

with Jo the total spin operator for the conduction system, referenced to the filled Fermi sea. 

That the Fourier components Q„ (n = — c» . . . oo) obey the level- /c Kac- Moody commutation relations with 
k — kg = 2 is easily verified by direct substitution of Eqs. (66-69): 

[Qn, Ql] - *e°''Qn+m + nSabS„. (71) 

Similarly, it is easily shown that the components Jn obey the level-fc Kac-Moody commutation relations with 
k = kj = 10: 

(72) 

Note that although the level of the axial-charge algebra, fc = 2 is equal to that for the two-channel impurity spin- 1/2 
Kondo problem |^, the level of the spin Kac-Moody algebra. A: = 10 is substantially higher. Thus, although the 
single-channel conduction spin-3/2, impurity spin 1/2 model and the two-channel impurity spin-1/2 model have the 
same conduction-band degeneracy, their low-energy spectra are considerably different. As pointed out by Sengupta 
and Kim |]l6| , the Kac-Moody algebra specified by Eq. ( fz^ ) is an example of the general result for conduction spin Sc 
which gives kj = 2SciSa + l)(2S'c + l)/3. 

The idea is now to write the Hamiltonian as a quadratic form in the axial charge (or isospin) currents and the spin 
currents. That this is possible is suggested by the Harmonic oscillator like structure of the commutation relation of, 
e.g., with Hcb- Specifically, 

[H,t,J?;^]^-VFm^J:i^. (73) 

By assuming such a separation into quadratic forms, directly computing the commutation and matching the nor- 
malization of the quadratic forms to Eq. (|7^), we can construct the Sugawara hamiltonian. It turns out that the 
normalizations depend only on the levels of the Kac-Moody algebras in the problem. 

The Kac-Moody levels determine the Sugawara form of the conduction Hamiltonian |pl| : 



He,- — 



^ ^ ~T • J— 71 ' Jn • + ^ ^ "7 • Q — n ' Qr 



12 ■ " " ^4 

n— — oo n— — oo 



(74) 



where : A : indicates normal-ordering of the operator A. 

We wish to diagonalize this Hamiltonian. To this end, we consider first the primitive states (or highest-weight 
states), which are eigenstates devoid of particle-hole excitations and satisfying the conditions |Q g < A:q/2 = 1 
and j < kj/2 = 5. We note that the combination of these constraints limits the primitive state spin values to 
j < 7/2. For n > 0, it follows from Eq. ( |6^ ) that the z component Qf^ annihilates particle-hole excitations, and 
hence Q^>ol0o) — for any primitive state \4>o). (For n > 0, by contrast, Qn\4>o > need not vanish, for Qn>o creates 
particle- hole excitations). Likewise, Qt>o\M = Qn>o\M = 0, since (Qn) creates (annihilates) electron pairs 
with one particle above the Fermi level and the other below. Finally, Jn>o annihlates particle-hole excitations and 
thus Jn>o\(t'o) ^ 0. It results that, if \q,j)o is the primitive state with axial charge q and spin j, then on the right-hand 



side of Eq. ( |74| ) only the n ~ term contributes to its energy Eo(q,j) 



Hcb\q.jh = ^ ( Y^^o + 4Q0 ) l9,J>o- (75) 



Since as Eqs. (|68[-[7lD show, Qo and Jo are the conduction-band axial charge and spin operators, respectively, the 

(76) 



primitive-state energies are given by 

"'Z('Z + 1) , j(.? + l) 



Eo{q,j) = —- 



12 
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To determine the quantum numbers q and j associated with each primitive state, we have to construct it from the 
single-particle levels of the conduction-band Hamiltonian. Particle-hole symmetry disposes these levels symmetrically 
with respect to the Fermi energy e p . For even number of levels (The number of levels depends on the number of lattice 
sites and on the boundary conditions; for periodic boundary condition, for instance, the number of levels equal the 
number of sites.), then, half of them lie above ep and half below, and the ground state of Hct will be nondegenerate. 
This spectrum corresponds to an odd numbered iteration of the NRG, and to antiperiodic free-fermion boundary 
conditions. For odd number of energy levels, in contrast, one level — which can accommodate up to four electrons — 
must lie at the Fermi energy, and the ground-state degeneracy is 2"* = 16. This corresponds to the even numbered 
NRG iterations and to free electrons with periodic boundary conditions. 

For an even number of levels, the nondegenerate ground state has, therefore, quantum numbers q = j = 0. The 
other primitive states are generated by filling the first level above ep with one or two electrons (or by vacating the first 
level below ep), corresponding to quantum numbers {q,j) equal to (1/2, 3/2), (1, 2), (0, 3), and another primary state 
corresponding to three elementary excitations, which is (1/2,7/2). All other excited states are not primary: since q 
increases as the energy levels are successively filled, the ones not containing elementary particle-hole excitations must 
have q > 1. The five primary multiplets appear in Table |. 

The 16-fold degenerate ground state comprises three multiplets: (i) the triplet {q,j) = (1,0), corresponding to a 
vacant (or fully occupied) level at ep, (ii) the octuplet {q,j) — (1/2,3/2), corresponding to a single electron (or three 
electrons) at ep, and (iii) the quintuplet {q,j) = (0,2), corresponding to a doubly occupied level at ep. Along with 
the states with two and three elementary excitations, {q,j) = (1,3) and (1/2,7/2), respectively, these are the only 
primitive states, since filling any level above the Fermi level (or vacating any level below ep) makes q > 1. They are 
listed in Table y. 

The eigenstates \(j)) of H^t containing one or more particle- hole excitations or with q > I are called descendant 
states; each such state can be obtained from a primitive state by successive applications of the raising operators Qn 
and J„ (n = 0,-1,-2,...). Disposed in order of increasing energy, a primitive state and its descendants form a 
conformal tower. 

Two eigenstates belonging to a same tower may have different quantum numbers q or j, and two energies on a 
same tower may have different degeneracies. The conformal towers are nonetheless important, for the commutation 
relations ([ti] ) and ( fz^ ) make the energy spacing in each tower uniform, the splitting AE — vptt/£ separating any two 
successive energies. It follows that, to calculate the low-energy spectrum of Hct, we need only construct the conformal 
towers associated with the lowest-lying primitive states. Eq. ( [76| ) gives the energy of the primitive state at the base 
of a tower, and the energies of the descendant states are given by 



q{q+l) , i(j + l) , 
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(77) 



where to is a positive integer. Notice that the q and j arc the quantum numbers of the primitive states at the base of 
the tower, in general different from the quantum numbers of the descendant states. The eigenvalues of the descendant 
states listed in Tables (|) and (||) were computed with Eq. (|77|). 



B. Interacting system 

We now consider the interaction term Hi. Expressed in terms of the spin Fourier components J„, the right-hand 
side of Eq. (|6^) takes a particularly simple form: 

oo 

Hi = ^pJSj- (78) 

n— — oo 

To take advantage of the similarity between this expression and the one defining the Sugawara Hamiltonian, Eq. ^ 
it is convenient to introduce the shifted Fourier components 

Jn = Jn+Sl. (79) 

Since the spins J an Sj have the same [SU(2)] symmetry, these shifted components follow the (level 10) commutation 
relations [Eq. (fz^)] obeyed by J„. Moreover, for pj = 1/6, Eqs. ( |7^ ) and (|7|) combine into a single expression for 
the model Hamiltonian (1) that is formally equivalent to the conduction Hamiltonian ([72[): 
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oo ^ oo - 

^1 = ^ E T^:^-«-Jn:+ E J-Q-n-Qn (80) 



71— — oo n— — C30 



Identifying this special coupling constant with the low-temperature fixed point of the model Hamiltonian, we reduce 
the computation of the low-energy sp ectru m to the diagonalization of a quadratic Hamiltonian. This diagonalization 



follows the steps outlined in Section III A, There, however, the quantum numbers q and j of the conformal towers 
were determined by the free-electron single-particle levels. Here, to determine the quantum numbers q' and j' of the 
interacting system, we refer to Affleck's and Ludwig's fusion rule Q. 

According to that rule, the embedding of the impurity spin Si = 1/2 in the conduction-band spin dictated by Eq. |7^ 
leaves the conformal tower axial charge unchanged, q' = q and yields a sequence of spins j' ranging from |j — 1/2| to 
the minimum oi j + 1/2 and kj — j — 1/2. Thus, for j < 9/2, i. e., for the spins we are interested in, the fusion rule 
is equivalent to the angular momentum addition law. 

It is therefore a simple matter to construct the conformal towers for the interacting system, the low-energy spectrum 
being given by Eg. ( p O |) . Results for the cases of a nondegencrate and of a degenerate free-electron ground state appear 
in Tables HI and p^V|, respectively. 

Finally, we briefly address an issue raised by Sengupta and Kim iQ, who suggest that for arbitrary Si,Sc, since 
the spin current Kac-Moody algebra is the same as for a multichannel model with kj = 2Sc{Sc + l)(2S'c -f l)/3 
channels that overcompensation will occur for impurity spins obeying Si < kj/2. This is based upon the result for the 
multichannel model which says overcompensation occurs for Si < k, k the number of channels. We disagree with this 
claim. To see why, we note that based upon the single fusion rule the maximal ground state spin which is obtained 
considering both non-degenerate and degenerate free spectra is 5/ (degenerate) and S* — S*/ (non-degenerate) where 



Sa 

E 

m>0 



_ c{ c + — ) even] (81) 



{Sc + 1/2) 



2 



[25c odd] . (82) 



In view of this, we believe the maximum overcompensated spin value is determined by the condition S* > Si, not 
kj/2 > Si. This checks with the generalization of the strong coupling picture presented in Fig. 2, since the largest 
conduction spin you can form from the single band in this picture is in fact S*, which restricts overcompensated 
impurity spins to Si < S* — 1/2. 



C. Boundary Operator Spectra and Physical Properties 

The conformal theory may also be used to generate the spectra of boundary operators about the non-trivial fixed 
point by evaluating the spectrum for a different set of boundary conditions, as noted by Affleck and Ludwig 
Specifically, one considers the situation where a Kondo impurity is placed at each end of the system. Then one applies 
a double fusion rule to identify the spectrum, because the impurity spin must be absorbed at each end. Denoting the 
impurities at either end by 5*/, S'j the double fusion rule that yields allowed S values of this double impurity system 
is, for a parent state with spin S' 

0<S<min{S' + Si + S'i,kj-S'-Si-S'i} . (83) 

The quantum numbers of the states with these boundary conditions then identify the quantum numbers of allowed 
boundary operators. These operators are so called because the critical behavior is limited solely to the boundary of 
the system. The normalized energy levels iE/vpT^ give the scaling dimensions Aq of the operators, which enter into 
time dependent correlation functions of the boundary operators through the long time expression 

< O(r)O(0) >^ . (84) 

We present the spectrum of primary field boundary operators in Table ^ as determined by the double fusion rule. 
From this table we can glean immediately a number of facts about the physical properties of the model: 
1) Low Temperature Behavior of the Specific Heat and Magnetic Susceptibility. The field with quantum numbers 
q = o, j = 1 is a triplet of primary spin field operators ^^-v-^, which has scaling index of = 1/6. In the manner of 
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Affleck and Ludwig we may identify the leading irrelevant operator about the fixed point Hamiltonian as ■ $. 
If we carry out second order perturbation theory to compute the impurity contribution to the specific heat coefficient 
C/T = 7 and magnetic susceptibility x precisely as did Affleck and Ludwig for the multichannel Kondo model 1^ , we 
will obtain 7,x ~ y^A^-i _ 2^-2/3^ This conclusion was also reached by Sengupta and Kim fl^ . 

2) Crossover in Applied Magnetic Field Application of a magnetic field is a relevant perturbation which, in analogy 
with the multichannel model, will drive the system to a Fermi liquid fixed point associated with a polarized scatterer 
p3[ . The scaling dimension of the magnetic field h is 1 — = 5/6 and the crossover scale Th = h'^'^ /t]I^ . To see 
this, consider adding the local Zeeman term 

Szeeman = h' / dr$^ (t) (85) 

Jo 

to the action. Upon rescaling time by A, the prefactor becomes h^A^^^° due to the singular behavior of $^ for long 
times. This must be a constant to preserve the form of the action. Hence (/i)^/'"'^^^* = (h^)^^^ ^ A and the physical 
properties will be universal functions in x = h^/^ /T . 

3) Relevance of Exchange Anisotropy The operator Q with quantum numbers g = 0, j = 2 is a local quadrupolar 
field. Because its scaling index is 1/2, this suggests that application of a uniaxial stress or electric field gradient 
conjugate to the quadrupolar field will be a relevant perturbation. This is even the case for Si = 1/2, 3/2 in contrast 
to the multichannel Kondo model. The reason is straightforward: exchange anisotropy in this model will split the 
conduction band degeneracy between jS"^! = 3/2, 1/2 states which is responsible for the nontrivial ground state. This 
is analogous to the lifting of channel degeneracy which is always relevant in the multi-channel Kondo model. In the 
event that we generalize from impurity spin 1/2 to impurity spin 1, the exchange anisotropy will also crystal field 
split the impurity spin and further alter the fixed point structure. 

4) Presence of Fermion Operators The operators with q = 1/2, j = 3/2 have appropriate quantum numbers and 
scaling indices to be identified with fermion fields. As in the two-channel model |^ there are two such fields here, and 
the interpretation of this is unclear, though it appears to be a generic result for overcompensated Kondo models with 
Si = 1/2. 



D. Comparison between NRG and Conformal Theory 



In the last columns of Tables III and IV are NRG eigenvalues normalized such that the first excited state energy 
relative to the ground state is forced to agree with the conformal theory finite size spectra. The scaled Energies are 
then multiplied by a factor of 24 to multiply out the lowest common denominator of the conformal theory spectra. 
NRG iterations for odd number of sites correspond to even number of levels and so should be compared with the 
conformal theory spectra for non-degenerate free electrons. NRG iterations with even numbers of sites correspond to 
odd number of levels and thus should be compared with the conformal theory spectra for degenerate free electrons. 

We can readily see from the tables that with few exceptions the agreement is quite satisfactory. When discrepancies 
arise, there are three sources, identified in Ref. (i) the logarithmic discretization of the NRG conduction band 

breaks conformal invariance, and will lead to detailed deviations between NRG and conformal theory spectra as 
the energy is raised; (ii) the conformal theory spectra themselves are generated at the fixed point, and corrections to 
scaling can arise as the energy is raised due to the irrelevant operators about the fixed point (the NRG takes these into 
account); (iii) the truncation of states in each q,j sector of the NRG will lead to systematic errors in the eigenvalues 
as one moves the excitation energy towards the cutoff value. 

It is satisfactory that in each sector (degenerate vs. non-degenerate free states) we not only have good agreement 
on the energies but also find all the states generated by the conformal theory up to the energies considered. Hence we 
have little doubt that the non-trivial fixed point identified by absorption of the impurity spin in the conformal theory 
is precisely the same as the fixed point identified by the NRG analysis. 



IV. DISCUSSION AND CONCLUSION. 



We have analyzed the Si — 1/2, Sc — 3/2 single band Kondo model with the Numerical Renormalization Group 
(NRG) and conformal field theory. Our work confirms that the model displays a non-trivial intermediate coupling fixed 
point which is unstable to the relevant perturbations of an applied magnetic field and spin exchange anisotropy. In 
particular, we have obtained detailed numerical agreement between the finite size spectra generated by the conformal 
theory with those of the NRG, both for degenerate and non-degenerate free electron spectra. From the conformal 
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theory, we infer that the impurity contributions to the specific heat coeflicient and magnetic susceptibihty will diverge 
as T~'^l'^ and that in the presence of a magnetic field the system will crossover to a Fermi liquid with a polarized 
scatterer, with crossover scale = h^^^ /t]1^ and that thermodynamic properties will be universal functions of 
X — h^/^/T. We have argued that when the model is generalized to arbitrary Si,Sc that non-trivial intermediate 
coupling fixed points corresponding to over-compensation will arise for Si < with 5* defined by Eqs. |l| and 
|8^ . We have also generated the operator primitive boundary operator spectrum of the model employing the double 
fusion rule. Our conformal theory results agree with those of Sengupta and Kim [T^ ] for the finite size spectra and 
low temperature thermodynamics despite the different choice of basis for the Sugawara hamiltonian in our work. 

We wish to close by discussing the experimental relevance of our work and further theoretical issues raised by it. 

Our work may be of relevance at intermediate temperature scales in two-level system Kondo impurity materials as 
an unstable fixed point. An experimental example of this may have been realized in metallic point contact devices 
p2| . Such fixed points were first studied by Zarand |l3| . Similarly, this fixed point may be relevant at intermediate 
temperature scales for dilute alloys of Ce"^+ ions, such as Lai_a;Ce2;Cu2Si2 |3^, such as argued by Kim |Q and Kim 
and Cox ^ . 

On the theoretical front, two key issues are raised. First, we are unable to compute the residual entropy and 
residual resistivity of this model using the same methods of Affleck and Ludwig because as far as we are aware no 
generalization of the Kac-Peterson |3^] and Verlinde |35j formulae exists for SU{2) representations other than spin 
1/2. The computation of the residual entropy and residual resistivity rests upon knowledge of the Kac-Peterson ]3^ ] 
and Verlinde formulae for the spin 1/2 representation of SU{2). Such a generalization would be of theoretical 
interest at least for working out the general properties of this model for arbitrary Si,Sc, though it is unlikely to be 
of great experimental relevance. 

Second, and of greater interest is the connection of this model to the Bethe-Ansatz solutions of the multi-channel 
Kondo model, which point was first made by Zarand ||l^ in his 1/fc expansion solution to the /c-channel two-level 
system Kondo model. In the approach of Tsvelik and Wicgman ||l^ to the fc-channel Kondo model, it is shown 
that a model with a single band of spin k/2 conduction electrons possessing an interaction to the impurity spin of 
a certain polynomial P{x) in x ^ Si ■ Sc{0) is in fact an integrable model. It is claimed that this model has the 
same Bethe-Ansatz spectrum in the compensated case as the corresponding fc-channel impurity Anderson model in 
the local moment limit where a Kondo description applies. It is then conjectured that in the ovcrcompensated case 
the single band spin k/2 model has the same spectrum as the fc-channel Kondo model. In the case S'/=l/2, P{Si ■ Sc) 
reduces to a simple Heisenberg exchange form. This claim has obvious physical appeal as illustrated in Fig. 1, since 
the strong coupling limit of the model indeed "fuses" the conduction electrons into a single large spin complex with 
magnitude fc/2. However, some doubt about the conjecture arises when one considers the non-interacting limit. In 
that case, there are 2fc free fermion branches to the fc channel model, whereas the single channel spin fc/2 band has 
fc -|- 1 free fermion branches. Namely, in the noninteracting limit, the number of degrees of freedom of the two models 
simply do not match. 

Andrei and Destrei 1 18 1 solved the Bethe-Ansatz directly on the multichannel model. In their approach, a dynamical 
fusing of conduction electrons into a spin fc/2 object describes the ground state. This again would appear to nicely 
match the strong coupling picture of Nozieres and Blandin |^ . The Bethe-Ansatz equations obtained by Andrei and 
Destrei are identical to those obtained by Tsvelik and Wiegman |17| on the basis of the conjecture described 
in the preceding paragraph. Indeed, both Bethe-Ansatz solutions yield thermodynamics and finite size spectra in 
excellent agreement with NRG and conformal theory results on the fc channel models, giving great confidence in the 
final Bethe-Ansatz equations, arrived at through completely different routes. An interpretation of the correspondence 
is that the Tsvelik and Wiegman [0 conjectured mapping amounts to an effort to match the dynamical fusion of 
Andrei and Destrei to a local Hamiltonian. This conjecture leads us to expect that the fc = 3 channel Si = 1/2 
model should have the same spectrum as the single band Sc = 3/2, Si = 1/2 model, but this is clearly not the 
case as shown here and by Sengupta and Kim |16|. We regard this discrepancy between the Tsvelik and Wiegman 
Bethe-Ansatz approach [|7| and NRG/conformal theory approaches presented here and in Sengupta and Kim |1^ as 
an interesting paradox to be resolved. 
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TABLE L Low-energy spectrum of the free conduction-band Hamiltonian for a nondegenerate ground state. Energies E 
(measured from the ground state) less than 2vftt/£ are shown. The fourth column distinguishes primary states (p) from 
descendant ones (d). 
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TABLE II. Low-energy spectrum of the free conduction-band Hamiltonian for a degenerate ground state. Energies E 
(measured from the ground state) less than VFiv/i are shown. The fourth column distinguishes primary states (p) from 
descendant ones (d). 
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TABLE III. Low-energy spectrum of the interacting model Hamiltonian for a nondegenerate conduction-band ground state 
(even iteration in NRG). Energies of the states derived from the ones in Table ^ by the fusion-rule procedure described in 
the text are shown. The fourth column distinguishes primary states (p) from descendant ones (d). The fifth column lists 
some NRG low energy spectra normalized such that the energy difference between the ground state and the first excited state 
agree with the conformal theory. Overall the agreement is excellent. The lack of perfect numerical agreement is due to the 
combination of the breaking of conformal symmetry by the NRG logarithmic discretization scheme and the truncation of states 
above a particular energy necessary to render the NRG block Hamiltonian matrices finite. Note the reversal of the pairing 
of integer(half-integer) q to integer(half-integer) j compared to Table ^ which indicates the non-Fermi liquid character of this 
spectrum in conjunction with the fractional energy spacings. 
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TABLE IV. Low-energy spectrum of the model Hamiltonian for a degenerate conduction-band ground state (odd iteration 
in NRG). Energies of the states derived from the ones in Table |l^ by the fusion-rule procedure described in the text are 
shown. The fourth column distinguishes primary states (p) from descendant ones (d). The fifth column lists some NRG low 
energy spectra normalized such that the energy difference between the ground state and the first excited state agree with the 
conformal theory. Overall the agreement is excellent. The lack of precise numerical agreement is due to a combination of the 
breaking of conformal symmetry by the NRG logarithmic discretization procedure and the truncation of energy levels to render 
the NRG block Hamiltonian matrices finite for diagonalization. Note the reversal of the binding of integer(half-integer)g to 
integer(half-integer) j in comparison with Table which, along with the fractional energy spacings, indicates the inapplicability 
of Fermi liquid theory to this spectrum. 
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TABLE V. Primary Field Boundary Operator Spectrum for the Sc = 3/2, Si = 1/2 model. We obtain this spectrum by 
applying the double fusion rule of AfHeck and Ludwig |^ to the spectrum of Table This rule assumes an impurity at 
each end of the finite size chain of length £. The scaling dimensions Ao are simply the scaled energies lE/vp-w for this double 
impurity boundary condition. The Q = 0,j = operator is simply the local charge of the site. The physical interpretation 
of the q = Q, j = 1 operator is the primary spin field at the impurity site. With the scaling dimension of 1/6, this indicates 
that the specific heat coefficient and susceptibility will diverge as T~^^^ as argued in the text. The q — 1/2, j = 3/2 field is a 
fermionic operator, which always must have scaling index 1/2 in any model. Fields with q = 1 include pair fields [q^ = ±1), 
which clearly do not provide singular pair field susceptibilities here. The q = 0, j = 2 field is a quadrupolar tensor, which given 
scaling index 1/2, indicates the relevance of exchange anisotropy for this model. See the text for further discussion. 
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FIG. 1. Strong coupling limit of the two-channel Si = 1/2 Kondo model At strong coupling (zero kinetic energy) 
in the two-channel spin 1/2 Kondo model, two units of conduction spin (labeled cl and c2) bind to form a net spin 1 complex 
that lines up antiparallel to the impurity spin. The resulting spin of the ground state is Stot = 1/2. The strong coupling 
fixed point is unstable to the introduction of the kinetic energy because antiferromagnetic superexchange will be generated to 
electrons off the impurity site as argued by Nozieres and Blandin Q of order t^/J, where t is the intersite hopping, which maps 
the effective model back to the weak coupling limit where the exchange interaction must grow. Since both weak coupling and 
strong coupling limits are unstable, a non-trivial fixed point must exist at intermediate coupling. 



FIG. 2. Strong coupling limit of the single channel Sc = 3/2, 5*/ = 1/2 model In the strong coupling limit of 
this model, two conduction electrons with Sz = 3/2, 1/2 will bind to form a spin 2 object that aligns antiparallel with the 
impurity spin. This strong coupling fixed point is unstable to the introduction of the kinetic energy which will generate an 
antiferromagnetic superexchange with the Stot = 3/2 bound object of order /J, which then maps it back to a weak coupling 
Kondo problem with growing exchange interaction. Since both strong- and weak-coupling limits are unstable, a non-trivial 
fixed point must exist at intermediate coupling. 



FIG. 3. RG flow diagram in the free electron limit for Sc = 3/2,5/ = 1/2 exchange interaction model. The RG 

fiow diagram is displayed for bare coupling g = 10"^^°. This bare coupling corrrespond to the free electron limit. The low lying 
exciattion energy spectrum is that of free conduction electrons. Note that low lying energy levels are evenly spaced. Labeling 
the first five lowest lying symmetry sectors are as follows throughout all the RG fiow diagrams presented here. (1) A*' = odd 
case, the ground energy symmetry sector is (0,1/2): the solid line; (1/2,1): the dashed line; (1/2,2): the dash-dotted line; 
(0,5/2): the dash-dot-dotted line; (1,3/2): the dash-dash-dotted line. (2) N — even case, the ground energy symmetry sector 
is (0,3/2): the solid line; (1/2, 1): the dashed line; (1,1/2): the dash-dotted line; (1/2,2): the dash-dot-dotted line; (0,5/2): 
the dash-dash-dotted line. 



FIG. 4. RG flow diagram displaying the crossover from the weak coupling to non-trivial finite coupling fixed 
point for Sc = 3/2, Si = 1/2 exchange interaction model. The RG fiow diagram is displayed for initial exchange coupling 
g = 0.01. With this weak coupling, we can observe the clear crossover from the weak coupling fixed point to new fixed point 
which is completely different from the former. Before the crossover, the low lying exciattion energy spectrum is close to that of 
free conduction electrons. After the crossover, the low energy excitation spectrum is described by non-Fermi liquid. See Figs. 
II, audi 



FIG. 5. RG fiow diagram with the initial weak coupling for Sc = 3/2, Si = 1/2 exchange interaction model. The 

RG fiow diagram is displayed for bare coupling g — 0.1. This bare coupling corresponds to the weak coupling regime. The low 
lying excitation spectrum is quite different from those of free electrons in Fig. 0. 



FIG. 6. RG fiow diagram with with the initial coupling close to the fixed point for Sc = 3/2, Si = 1/2 exchange 
interaction model. The RG flow diagram is displayed for bare coupling g — 0.2. This bare coupling is close to the fixed 
point. Note that the third-order scaling analysis gives g* = 0.2 at the fixed point. 
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FIG. 7. RG flow diagram with the initial strong coupling for Sc = 3/2, Sj = 1/2 exchange interaction model. 

The RG flow diagram is displayed for bare coupling g = 0.5. This bare coupling corresponds to the strong coupling regime. 
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